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Abstract 

Inclusive semileptonic decays are discussed in the framework of heavy quark 
effective field theory by employing the short distance expansion in the effective 
theory. The lowest order term turns out to be the parton model; the higher 
order terms may be regarded as correction terms to the parton model result. 
The first nonvanishing corrections to the parton model result are given and 
the lepton energy spectrum of inclusive semileptonic decays of heavy mesons 
is calculated. 



1 Introduction 



Heavy quark effective field theory |], || has turned out to be a useful tool for the 
description of heavy quark systems. Due to additional symmetries Q emerging in 
the heavy quark limit several predictions for the decays of heavy hadrons may be 
obtained which are completely model independent. In particular, the form factors 
of weak decays involving a heavy to heavy transition are severely constrained by 
heavy quark symmetry and some predictions may be obtained even for heavy to 
light deacays [f|. 

The heavy quark expansion has also been applied to inclusive decays. The first 
discussion of inclusive semileptonic decays in this framework has been given by Chay, 
Georgi and Grinstein in Q. Along similar lines proceeds the work of Bigi, Shifman, 
Uraltsev and Vainstein || and of Blok, Koyrakh, Shifman and Vainshtein || , where 
inclusive nonleptonic as well as semileptonic decays have been considered. 

Usually the inclusive decay rates of heavy hadrons and their spectra are ap- 
proximated by the decay of the heavy quarks in a parton model approximation. 
However, this is a model approach and an estimate of its error is difficult. This 
in is contrast to the methods described in [|], [| and also to the present approch 
which is a controllable approximation to QCD corresponding to expansions in small 
parameters. 

The method proposed in ]4], [5| to deal with inclusive processes is to use an oper- 
ator product expansion, which is justified for heavy hadrons due to its large mass. 
The lowest dimension operator is a dimension three operator; its matrix elements 
are normalized due to heavy quark symmetries and yield the parton model result. 
Corrections enter through operators of dimension five; their matrix elements are 
parameters to be taken from measurements. These corrections have been calculated 
for nonleptonic and for semileptonic decays of mesons in || and for the differential 
distributions in semileptonic decays in ||; semileptonic decays of baryons have been 
considered in [^|. 

In the present paper we shall consider only semileptonic decays and formulate a 
slightly different approach as the one proposed in [pi H]. After switching to heavy 
quark effective theory the space time dependence of the heavy quark fields in the 
hadronic tensor is only due to the residual momentum k of the heavy quark Q 
inside the heavy hadron H, i.e. k = Ph — vhqv. In semileptonic decays there is 
in addition to the scale set by the heavy mass tuq also the momentum transfer to 
the leptons q 2 . After switching to the effective theory the relevant variable becomes 
Q 2 = (q — itiqv) 2 , which defines in addition to the mass ttiq of the heavy quark a 
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second scale in the problem. 

In large portions of the phase space Q 2 is of the order of mq and thus a simul- 
taneous expansion in both 1/rriQ and 1/Q 2 is appropriate. This is the approach 
chosen in [|5], However, one may also consider the region in phase space, where 

™ 2 Q > (? - m Q v ) 2 > ^QCD 

and thus an expansion only in powers of 1/Q 2 is necessary while only the leading 
term of the I/vtiq is kept. 

In the present paper we shall consider the latter approach and compare to the 
one of |5], |J. The paper is organized as follows. In the next section we fix our 
notation by giving the basic formulae for the inclusive decay rates. In section 3 we 
set up the operator product expansion in the effective theory and discuss the lowest 
order term and the first nonvanishing correction. The setup given here may be used 
for baryonic as well as for mesonic decays and in section 4 we give a parametrization 
for the hadronic matrix elements for meson- and baryon inclusive decays and discuss 
short distance QCD corrections. Finally, in section 5 we discuss the lepton energy 
spectrum in inclusive semileptonic decays in some detail. 



2 Inclusive Decay: Matrix Elements and Rates 

In order to set up the notation we shall briefly recall the basic formulae for the 
description of inclusive semileptonic decays. The inclusive semileptonic decay of a 
heavy hadron H (either a heavy meson or a heavy baryon) into some hadronic state 
X and two leptons 

H(P H = m H v) -> X{P X ) + £(k) + Vi{k') (1) 
is mediated by an effective Hamiltonian of the form 

Heff = ^V Qq (Q^(l - 7b)ff)(^(l - Tb)"i) (2) 

where Q is the heavy quark contained in the in the heavy hadron f/Q. The differential 
decay rate may be written in terms of the hadronic tensor W^ v {q, v) 

W^(q,v) = f x {2-Kf5\P H -q-P x ) (3) 

<H(v)\(Q 1 , (l - 75 )g)|X >< X\(q ltl {l - j 5 )Q)\H{v) > . 

1 I thank Mark Wise for a clarifying dicussion on this point. 

2 We shall assume that the flavor of q may be determined from the hadronic state X. 
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We shall consider only the case where the spins of the final state leptons are not 
measured; hence the leptonic tensor is given by 



(4) 



The differential decay rate in the rest frame of the heavy hadron is then 



dT 




V Qq \ 2 W fMV (k + k',v)A^d(PS) 



(5) 



where d(PS) is the phase space integral 



d(PS) = / dkdk <5(observables) dk 




d 4 k 



8(k 2 -m 2 )Q(k ). 



(6) 



The delta function <5(observables) projects out the observables to be considered; 
for the case of the energy spectrum of the charged lepton in the rest frame of the 
decaying heavy meson it is ^(observables) = 5{E — vk'), since k! is the momentum 
of the charged lepton. Furthermore, rriLep is the lepton mass which we shall neglect 
in the following. 

The hadronic tensor may be rewritten using standard techniques 



where we have explicitly spelled out the x dependence of the field operators. 

The tensor W^ u may be decomposed into scalar form factors which then depend 
only on two invariants; these may be cenveniently chosen to be q 2 and vq. As 
discussed in some detail in j| , one may relate these form factors to the disconitinuity 
of the ones defining the time ordered product 
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dV x < H(v)\(Q(x)^(l - 7fc)g(aO)(ff(0)7*(l - k)Q(0))\H(v) > 



(7) 




(8) 



<H(v)\T (Q(x)^(l- l5 )q{x))(q(0)l»(l-l5)Q(0)) \H(v) > 



For a fixed value of q 2 there are two cuts of (H) in the complex vq plane. A physical 
cut extents along the real axis for 




(9) 
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where /x is the mass of the lightest hadronic final state. In addition, there is a cut 
extenting along the real axis for 

vq > -^—{{2m H + fi) 2 - m 2 H - q 2 ) (10) 

The problematic region in phase space is the one where the invariant mass of the 
hadronic final state approaches \x\ here resonances will dominate and perturbative 
QCD will receive large corrections. However, even in this region the predictions 
of perturbative QCD will be recovered by a suitable smearing, i.e. by calculating 
a smooth average along some portion of the cut. This has been discussed in some 
detail for the case of e + e~ — ► hadrons in || and applied to the present case in ||. 
In particular, one may expect a suitable smearing, if the energy spectrum of the 
charged lepton is calculated. 

In order to calculate such an average one has to perform a contour integration of 
in the complex vq plane with a suitable weight function. The contour encircles 
the physical cut in such a way that the discontinuity T^ v (i.e. the weighted hadronic 
tensor W^) of the cut is picked up. As has been argued in [|J one may deform 
the integration contours in the complex vq plane in such a way that one always 
stays away from the dangerous region by a distance large compared to Aq C d- This 
may be achived provided that ji ^> Aq C d because for \i — > the cut ([]) and fllUp 
pinch the integration contour for q 2 = q^^ = {m H — /i) 2 . From these arguments 
one expects that perturbative QCD may be applied after suitable smearing and 
for jj, ^> Aqcd which means that we may only study inclusive semileptonic b —>■ c 
transitions in perturbation theory, while b —>■ u transitions are likely to receive large, 
nonperturbative corrections, at least in the region, where q 2 is close to its maximal 
value. 

Keeping this in mind one may relate the perturbatively calculated time ordered 
product ([§D to the hadronic tensor using the relation 

W^(q,v) = 2lmT fMU (11) 

The matrix element of the time ordered product in (|j) may be expansed in a 
sum of local operators by performing a short distance expansion. However, there are 
still two large momenta involved, the momentum of the heavy meson which scales 
with itlq and the momentum transfer to the leptons q 2 . In order to disentangle long 
and short distance contributions it is convenient to switch to the effective theory. 
In this process the heavy hadron momentum is rescaled and only the small residual 
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momentum remains. As will become clear below, the relevant momentum transfer 
variable will be Q = q — tuqv which we assume to be the only large scale in the 
problem. Once we are working in the effective theory, an operator product expansion 
will be an expansion in inverse powers of Q 2 . 

3 Operator Product Expansion in the Effective 
Field Theory 

Switching to the effective field theory allows to scale out the large parameter m#. 
The relevant momentum variable characterizing the initial state is then a small 
momentum, namely the residual momentum k = Ph — viiqv of the heavy quark. 
Matching the effective and the full theory to lowest order in the 1/rriQ and a s (mg) 
expansion one has to replace 

Q(x) - e- im ^-^h v {x) (12) 

where the x dependence of the static heavy quark field h v is due to the residual 
momentum k. 

Injecting this into the hadronic tensor one obtains 

W, u (q,v) = jd A xe ix ^- m ^ (13) 
< H(v)\(h v (x) lfl (l - l5 )q(x))(q(0)l^ ~ ls)h v (0))\H(v) > 

An operator product expansion (i.e. a short distance expansion) is justified, if 
the momentum transfer variable Q = q — rrigv is large compared to the residual 
momentum k. In addition, switching to heavy quark effective theory implies an 
expansion in powers of I/tuq. In contrast to [^|, we shall consider the portion 
of phase space in which Q 2 <C mn and thus the lowest order term in the I/tuq 
expansion is sufficient. Hence we may use ( |TB| ) as the starting point. 

Using momentum conservation Ph = rriQV + k = q + Px we find that an operator 
product expansion is possible for P\ ^> m 2 x and that this expansion is a power series 
expansion in powers of Aqcd/^x, up to logarithmic terms which may be resummed 
by renormalization group techniques. Note that, if mx is much smaller than m#, the 
use of the effective theory to lowest order is justified, since higher order corrections 
to the effective theory will be of the order Aq C D/ m H an d thus will be much smaller. 
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However, as argued above, we shall only consider b — > c inclusive transitions and 
thus mx Aqcd. For this case we may use the framework of perturbative QCD 
to calculate the coefficients of the operator product expansion. 

Thus we are lead to consider the time ordered product of the two currents 

T, v (q,v) = J d 4 xe* xQ T [(h v (x)^(l - ^ 5 )q(x))(q(0)^ - 7s)^(0))] (14) 

where Q = q—rriQV is the momentum transfer variable defining the large momentum. 
We shall define the operator product expansion for the projection of T^ v needed to 
obtain the inclusive semileptonic rate 

T(k,k',v)=T flu (k + k',v)A^ (15) 

Before we write down the general form of the short distance expansion in the effective 
theory we shall recall some basic formulae concerning the spin structure in the 
effective theory. The sixteen covariants built from Dirac matrices reduce to only 
four when projected into a definite velocity sector. We shall use as a basis the unit 
matrix and the three spin vectors s M with v ■ s — which correspond in the rest 
frame v = (1, 0, 0, 0) to the three Pauli matrices. The translation table of the sixteen 
invarants is given by 

1 — P+ = \{l + f) 1^ P+1»P+ = v„P + 

7^75 — ► P+1^P+ = s M 75 — ► ^+75^+ = 

The spin vectors satisfy the relation 

s^s v = (-g^ + u„u„)l + ie aitv pv° 's 13 (16) 

which is a generalization of the well known relation between the Pauli matrices and 
which may be used to reduce products of spin vectors. 

At tree level one may obtain the operator product expansion up to operators of 
dimension five by considering the two diagrams gepicted in fig.l. The result is 



T(k,k',v) = ( — 1 \u{k'Q)h4{l- l5 )h v (17) 
Q z — mf + te J 



1 

<2 2 — mf + ie 



+ (o2 TTT, ) 64^(1^)^(1 - lb )h v 
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I 128(k'Q)h v (tQD)Ul- l5 )h v 
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\2 



64(^2)^(^)^(1-75)^ 
Mh v [{iQD),{ik'D)\lJ-{l- lb )h v 
me Xa p p Q x k' p h v (iD a )(iD p )^l - l5 )h v 
128h v (ik'D)(iQD)^(l - -f 5 )h v 

+ (752 Vr^) 256(^2)^(^^)^(1-75)^ 

y - + ie J 

+ Operators of Dimension six or higher 
We note that the spin structure of all the terms is of the form 

KlfiO- ~ lb)K = v^h v h v - h v s^h v (18) 

reflecting the left handedness of the current. 

The differential rate may now be calculated from the imaginary part of the 
forward matrix element of the operator T 

dT = — -\V Qq \ 2 {2 Im < H(v)\T(k,k',v)\H(v) >)d(PS) (19) 
Amu 

The relevant hadronic matrix elements will be discussed in the next section. 



4 Hadronic Matrix Elements and Short Distance 
QCD Corrections 

We shall first consider the decay of a pseudoscalar meson. The only nonvanishing 
matrix element of the dimension three operators is 

< H(v)\h v h v \H(v) >= 2m H . (20) 
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which is normalized as given above due to heavy quark symmetries. Inserting this 
into ( |19| ) reproduces the parton model result. 

All the matrix elements of the dimension four operators vanish. The general 
form of the matrix element not involving the spin vector is 

< H(v)\h v D a h v \H(v) >= B v a (21) 

and by contraction with v and using the equations of motion we have B = 0. The 
corresponding matrix element involving the spin vector vanishes due to parity. 

The matrix elements of the dimension five operators may be expressed in terms 
of two independent parameters. Due to the equations of motion we have 

< H{v)\h v {iD a ){iD p )h v \H{v) >= -ImnB^g* 13 - v a v p ) (22) 

where B^ is determined by the trace 

< H(v)\h v (iD) 2 h v \H(v) >= -§m n B x . (23) 

This trace is related to the expectation value of the residual momentum squared 
k 2 of the heavy quark in the heavy meson (< k 2 >~ 3Bi) which is expected to 
be of the order of Aq CD . Unfortunately one may not extract a precise value from 
experimental data and one has to rely on theoretical estimates. An estimate based 
on QCD sum rules |10[ yields reltively large value of B% ~ 0.3 GeV 2 . We shall use 
a default value of B\ ~ 0.2 GeV 2 and discuss the dependence of the result on this 
parameter. 

The other independent matrix elelement is parametrized as 

< H{v)\h v {iD a ){iD p )s p h v \H{v) >= 2m H B 2 e Xa/3p v x (24) 

where the constant B 2 is given in terms of the scalar matrix element 

ie XaPp v x < H(v)\h v (iD a )(iDfs)s p h v \H(v) >= ig < H [v)\h v a ^G pu h v \H [v) > (25) 

where g is the strong coupling constant. This matrix element is related to the mass 
splitting between the pseudoscalar meson and its spin symmetry partner vector 
meson 

ig < Hiv^Ka^G^h^Hiv) >= -l2m H (m 2 B * - m 2 B ) (26) 

3 This parameter is - up to a factor 3 - the parameter Ai as defined in the paper by Falk and 
Neubert g 
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from which we obtain 

B 2 = ~(m 2 H * - m 2 H ) ~ 0.07 GeV 2 (27) 
8 

Thus lowest order corrections in the mesonic sector are given in terms of two pa- 
rameters, one of which is related to an observable quantity. 

The decay of a heavy ground state baryon A& has been considered in Q in 
the approach proposed in ||. For this decay we have two matrix elements of the 
dimension three operators of which the normalization is known due to heavy quark 
symmetries. We have 

< A H (v,s)\h v h v \A H (v,s) > = 2m H (28) 
< A H (v, s)|/i tJ s At ^|A ff (t),s) > = 2m H s fl (29) 

where s M is the direction of the spin of the Ah baryon in its rest frame. 

As in the case of the mesons the matrix elements of the dimension four operators 
vanish due to parity and the equations of motion. 

Finally, the matrix elements of the dimension five operators are given in terms 
of only one parameter 

< A H (v, s^h^iDaj^D^h^A^^) > = -2m H Ax(g a p - v a vp) (30) 
< AHfas^hvtfDa^iDpjSphvlAfas) > = -2m H Axs^g al3 - v a vp) (31) 

since the heavy quark spin symmetry relates these two matrix elements. 

The parameters A%, B\ and B 2 parametrize the long distance QCD contributions 
to the inclusive semileptonic decay rates. Since we are dealing with an operator 
product expansion the short distance QCD contributions may be calculated as the 
scale dependence of the Wilson coefficients, the tree level value of which may be 
read off from fll7l) . 

Short distance corrections may be considered using the usual machinery of the 
renormalization group. The calculation of the short distance QCD corrections may 
be found in the literature. The anomalous dimensions of the scalar operators defining 



the constants A 1 , Bx and B 2 has been performed in ||11||. 
Due to reparametrization invariance |T2|] the operator Ox 

O x = h v (iD) 2 h v (32) 

has vanishing anomalous dimension and thus Bx as well as Ax are scale invariant 
quantities. The second operator 2 

2 = igh v a^G» v h v (33) 
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has nonvanishing anomalous dimension but does not mix with any other operator. 
From the renormalization group equation with the anomalous dimension from [] I J 



= -^*s(n)B 2 (n) (34) 

one obtains 

/ a (a) \ 9 /( 33 - 2 ™/) 
%)=FPt B 2 (m H ). (35) 

\a s {m H ) J 

with rif = 4 for B hadron decays. This scale dependence introduces an Q 2 depen- 
dence, since this is the relevant momentum transfer variable in the present case. 
The effect of these short distance corrections will be discussed for the example of 
the lepton energy spectrum in inclusive semileptonic decays of heavy mesons. 



5 An Application: Energy Spectrum in Heavy 
Meson Decays 

As an obvious application we shall consider the lepton energy spectrum of the 
charged lepton in an inclusive semileptonic decay of heavy mesons. To be spe- 
cific we shall consider the decay B —>■ X c e~v e Here X is a hadronic final state 
containing charm. 

The differential decay rate is given by 

fir c 2 r 

!^ = JfJL\V Q J 2 dkdk5(Q 2 -m 2 )5(E-vk')(2Im < H(v)\T(k,k',v)\H(v) >) 
dE AiriH J 

(36) 

The calculational details of the phase space integration are not very iluminating 
and are given in an appendix. However, unlike in the usual case of phase space 
integrals we have due to the relations 

2ttnf -j \ 2 , A = -(2 7 r)5'(Q 2 -m 2 ) (37) 



Q 2 — mf + ie 



2 



Mm U- m? + J = Wm-ml) (38) 



1 \ 3 



2ilm — ^— = -n5"(Q 2 -m 2 ) (39) 

\Q 2 -mf + it 
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also integrations involving the first and second deriavtive of the delta function 8' 
and 8". Flipping the derivative yields derivatives of the integrands. Since the QCD 
short distance corrections imply a dependence of B 2 on the kinematic variables one 
has to differentiate B 2 as well. This yields an additive correction term which is 
proportional to a s taken at the scale Q 2 ; this scale is fixed to Q 2 = mf due to the 
delta functions showing that the present approach is indeed limited to mi ^> Aq C d- 
The results for the lepton spectrum is displayed in fig. 2 and 3. Fig.2 shows the 
charged lepton energy spectrum for the inclusive semileptonic decay of a B° or a 
B~ into charmed final states. The parameters used are m# = = 5.28 GeV, 
mi =m c = mD = 1.86 GeV, \V cb \ = 0.043 and B 1 = 0.2 GeV 2 . The solid line (a) is 
the full result, including long- and short distance corrections. This is compared to 
the spectrum including the corrections only at tree level, i.e. without the QCD short 
distance contribution (line (b)). Finally, curve (c) is the leading order result which 
is the spectrum obtained from the parton model. The spectrum obtained this way 
is slightly shifted to lower energies and the short distance corrections yield a small 
enhancement at low energies. Qualitatively this result agrees with the one obtained 
in§. 

Finally, we shall investigate the dependence on the parameter B\ by considering 
the total rates. In fig. 3 we plot the quantity 

-f" corri^B -V c 6 Z^ e ) _^ 

^partoni^B > X C G Ve) 

in percent as a function of the value of B\. We chose the same parameters as for fig.2, 
except that Bi is varied in the range 0.04 < B 2 < 0.4. Curve (a) is the full result, 
including long- as well as short distance contributions, curve (b) conbtains the tree 
level corrections only. The corrections reduce the total rate (compared to the parton 
model rate) for values of B\ above 0.1. This seems to favor large values of Bi, since 
a reduction of the semileptonic branching fraction is phenomenologically welcome. 
Furthermore, one could in turn try to obtain matrix element corresponding to B\ 
from a preceise measurement of the inclusive lepton spectrum. 

6 Conclusions 

Inclusive semileptonic decays of heavy hadrons may be described by employing an 
operator product expansion for the operators appearing in the hadronic tensor. The 
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expansion is performed after switching to heavy quark effective field theory where 
the large momentum of the heavy hadron is scaled out. The relevant momentum 
transfer variable is Q = q—rriQV where q is the momentum transferred to the leptons. 

In the present paper the region of phase space has been considered in which 
Q 2 -C ttlq and heavy quark effective field theory to leading order in the I/tuq 
expansion is sufficient. The approach of || and || is in fact valid in a different 
portion of phase space, since || || assume that Q 2 ~ rnq and thus a simultaneous 
expansion in 1/ vtlq and 1/ Q 2 is necessary. 

The lepton spectrum as obtained in ||, || contains delta functions and its deriva- 
tives which reflects the failure of the expansion at the edge of phase space E ~ E max . 
In contrast to that the approach applied here yields a completely smooth lepton en- 
ergy spectrum without any delta function singularities. 

Of course, the limitations of the procedure proposed here (namely Q 2 -C ttiq) 
do show up at some point. If one studies for example the rate differential in the 
momentum transfer q 2 = (k + k') 2 and the total leptonic energy vq = vk + v k! then 
one finds also delta functions and its derivatives entering the differential rate using 
the method chosen here. However, after an integration over the total leptonic energy 
one obtains a smooth q 2 spectrum and the only reminder for the breakdown of the 
Q 2 expansion is a behaviour of the correction term of the type 

dr i 

T~2 ~ / ( 4 °) 

q V \ m Q - m l) 2 - <? 2 

at the endpoint of the leptonic q 2 spectrum. This is an integrable singularity and 
hence a finite total rate may be obtained as expected. 

In conclusion, depending on the region of phase space considered one has to 
choose between the method proposed in |5], |J and the one presented here. To ob- 
tain the lepton spectrum one has to integrate over the phase space and phase space 
regions contribute where the approximations used become invalid. However, both 
methods yield finite total rates which means that the singularities in the lepton spec- 
trum are intagrable and one may still expect a reasonable lepton energy spectrum 
from both methods; this expectation is confirmed by the qualitative agreement of 
the lepton energy spectra obtained from both approaches. 
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A Some Details of the Calculation of the Lepton 
Energy Spectrum 

In this appendix we give some of the technical details of the calculation for the 
energy spectrum of the charged lepton. 

In order to perform the phase space integrations in ( j36|) it is convenient to define 
three integrals of some function F of the three independent scalar products v k, vk! 
and 2kk' 

F = F(vk } vk' } 2kk'). 

We define the integral 

I [F] = J dkdk'5{Q 2 -rnf)6(E - vk')F(vk,vk',2kk') (41) 

P 



— J dkF(k, E, 2m H (E - E max + k')) 



32vr 3 

a 

where the limits of the integration are given by 

a = m H E ™~ 2 E E P = E max - E (42) 

By definition E is the energy of the charged lepton in the rest frame of the decaying 
heavy meson and E max is its maximal energy 

F m ax = 7T—( m H - m2 l) ( 43 ) 

The corrections contain derivatives of delta functions and it is convenient to define 
two more integrals 

h[F] = J dkdk'5'{Q 2 -m*)5{E - vk')F(vk,vk',2kk') (44) 

I 2 [F] = JdkIk5"(Q 2 -mf)5{E -vk')F{vk,vk',2kk') (45) 
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where 5' and 8" denotes the first and second derivative of the delta function. These 
two integrals may be expressed in terms of the integral Iq 

= l(F<J>,E,0)_F( a ,E,4Ec)\_ 
[ J 64vr 3 \ m H m H -2E J 01 J v ; 

UF] = ^-/(2 mH (E m - g))F(0, E, 0) —^ _ m + J„n (47) 

where F" and F" denotes the first and second derivative of F with respect to the 
last argument 2kk'. 

Inserting the operator product expansion fllTD one finds that the piece containing 
the dimension three operators reproduces the parton model result 

— = -32G 2 F \V Qq \ 2 I [(vk)(kk' - m H vk')\ (48) 

At tree level the correction terms originating from the dimension five operators 
are given by 

= 32G 2 F \V Qq \ 2 {B 1 I 1 [5(vk)(k'Q)-2(vk)(vk r )(vQ)] (49) 



dE 



+ B 2 h[2{vk'){kQ) - 2{kk'){vQ)}} 

QAG 2 F \V Qq \ 2 Bil 2 [{k' Q)(vk)((vQ) 2 - Q 2 )} 



The QCD short distance contributions may be included by putting in the scale 
dependence of B 2 as discussed in the last section. At the one loop level one finds 
that aside from the fact that the value of B 2 is changed additional terms originating 
from the derivative in ( |4*6D acting on B 2 occur, since now B 2 depends on 2kk! . The 
derivative of B 2 is obtained from the renormalization group equation fl34|). Inserting 
this we obtain an additional term due to the short distance contribution to B 2 

AsD JE = - 32G Fl^| 2 /o[2(^')(fcQ) - 2(kk , )(vQ)}^a s (m 2 )B 2 (m 2 ) (50) 

Here and in the other terms where it occurs B 2 has to be scaled down to Q 2 which 
is due to the delta functions of phase space fixed to the value mf. 
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Figure Captions 



Fig.l Feynman Diagrams to be evaluated to obtain the Wilson coefficients of the 
operators up to dimension five. The blob denotes the left handed current and 
the double line denotes the decaying heavy quark. 

Fig. 2 The lepton energy spectrum dT/dE x 10 15 of the inclusive semileptonic decay 
B X c e~v e . The parameters used are m& = mg = 5.28 GeV, m c = mo = 
1.86 GeV, \V cb \ = 0.043 and B x = 0.2 GeV 2 . The curves are: (a) Full result 
including long- and short distance contributions, (b) Result without short 
distance QCD corrections, (c) leading order result, i.e. the parton model 
result. 

Fig. 3 The ratio of total rates (r corr /rp arton ) — 1 in percent for the inclusive decay 
B° — > X c e~v e function of the parameter B x . The curves correspond to: (a) 
full result, (b) corrections only at tree level. 
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Fig.l: Feynman Diagrams 




Fig. 2: Lepton Energy Spectrum in Inclusive b --> c Decays 
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Fig. 3: Ratio of Total Rates 
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